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Abstract. We prove that when subjected to periodic forcing of the form p M , P)W (t) = 
p(ph(x, y) + sin(u;i)), certain second order systems of differential equations with dissipative 
homoclinic loops admit strange attractors with SRB measures for a set of forcing param- 
eters (p,p,uj) of positive measure. Our proof applies the recent theory of rank one maps, 
developed by Wang and Young [30, 34] based on the analysis of strongly dissipative Henon 
maps by Benedicks and Carleson [4, 5]. 



1. Introduction 

In this paper we establish connections between a recent dynamics theory, namely the 
theory of rank one maps, and a classical dynamical scenario, namely periodic perturbations 
of homoclinic solutions. We prove that when subjected to periodic forcing of the form 
P^ P ,uj(f) = n(ph(x, y) + smut), certain second order equations with a dissipative homoclinic 
saddle admit strange attractors with SRB measures for a positive measure set of forcing 
parameters (fi,p,uj). 

A. The theory of rank one maps. The theory of rank one maps, systematically developed 
by Wang and Young [30, 34], concerns the dynamics of maps with some instability in one 
direction of the phase space and strong contraction in all other directions of the phase 
space. This theory originates from the work of Jackboson [13] on the quadratic family 
f a ( x ) = 1 — o,x 2 and the tour de force analysis of strongly dissipative Henon maps by 
Benedicks and Carelson [5]. 

The theory of ID maps with critical points has progressed dramatically over the last 30 
years [19, 13, 8, 4, 28]. The breakthrough from ID maps to 2D maps is due to Benedicks 
and Carleson [4, 5]. Based on [5], SRB measures were constructed for the first time in [6] for 
a (genuinely) nonuniformly hyperbolic attractor. The results in [5] were generalized in [20] 
to small perturbations of Henon maps. These papers form the core material referred to in 
the second box below. 



Theory of 
ID maps 



Henon maps 
& perturbations 



Rank one 
attractors 



All of the results in the second box depend on the formula of the Henon maps. In going 
from the second box to the third box, the authors of [30] and [34] have aimed at developing 
a comprehensive chaos theory for a nonuniformly hyperbolic setting that is flexible enough 
to be applicable to concrete systems of differential equations. 
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The theory of rank one maps has been applied to various systems of ordinary differential 
equations [31, 32, 16, 11, 21]. The most siginificant application thus far has been the analysis 
of periodically- kicked limit cycles and Hopf bifurcations [31, 32]. In these cases, periodic 
kicks of the limit cycle and separated by long periods of relaxation to the limit cycle. If 
the contraction to the limit cycle is weak and the shear is strong, then admissible families 
of rank one maps are produced. The analysis of periodically-kicked Hopf limit cycles has 
been extended to the setting of parabolic partial differential equations [17]. These studies 
illustrate that the theory of rank one maps can be used to rigorously prove the existence of 
strange attractors with SRB measures for physically meaningful differential equations. 



B. Periodically-perturbed homoclinic solutions. Periodically-forced second order systems, 
such as the periodically-perturbed nonlinear pendulum, Duffing's equation, and van der 
Pol's equation have been studied extensively in the past [3, 9, 10, 14, 15, 29]. When a 
given second order equation with a homoclinic saddle is periodically perturbed, the stable 
and unstable manifolds of the perturbed saddle intersect transversely within a certain range 
of forcing parameters, generating homoclinic tangles and chaotic dynamics [18, 22, 23, 24, 
27]. Homoclinic tangles were first observed by H. Poincare [22, 23, 24]. There also exist 
parameters for which the stable and unstable manifolds of the perturbed saddle are pulled 
apart. For these two cases, Fig. 1 schematically illustrates the time-T maps for the perturbed 
equations, where T is the period of the perturbation. The first picture leads to exceedingly 
messy dynamics and the second appears simple. 




Fig. 1 (a) Homoclinic intersections and (b) separated invariant manifolds. 



In this paper we study periodically-perturbed second order equations but we follow a 
new route. Instead of looking at the time-T maps, we extend the phase space to three 
dimensions and we explicitly compute return maps induced by the perturbed equations in 
a neighborhood of the extended homoclinic solution. To be more precise, we use variables 
(x, y) to represent the phase space of the unperturbed equation and we let (x, y) = (0, 0) be 
the saddle fixed point. Write the homoclinic solution for (x, y) = (0, 0) as t. We construct 
a small neighborhood of I by taking the union of a small neighborhood U £ of (0, 0) and a 
small neighborhood D around £ outside of Ui £ . See Fig. 2. Let a ± G U e H D be the two line 
segments depicted in Fig. 2, both of which are perpendicular to the homoclinic solution. We 
use the angular variable 9 e S 1 to represent the time. 
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Fig. 2 U £ , D and a ± . 
In the extended phase space (x,y,9) we define 

U £ = U £ x S\ D = D x S 1 

and we let 

E± = a ± x S 1 . 

Let N : S + — > E~ be the map induced by the solutions in U £ and let M : £~ — > S + be the 
map induced by the solutions in D. See Fig. 3. We first compute M and USf separately. We 
then compose N and M to obtain an explicit formula for the return map No M : S _ — > 
We show that for a large open set of forcing parameters, these return maps naturally fall 
into the category of the rank one maps studied in [30] and [34]. 




Fig. 3 Construction of return maps. 

C. A brief summary of results. Autonomous second order systems with a dissipative 
homoclinic saddle are subjected to periodic forcing of the form p^ PtLU (t) = n(ph(x,y) + 
sin cut), where /i, p, and ui are forcing parameters. We prove that if the saddle is dissipative 
and nonresonant (see (HI) in Section 2) and if the unperturbed equation satisfies certain 
nondegeneracy conditions (see (H2) in Section 2), then there exists an interval [pi, P2} such 
that for p G [pi , po\ , the family of return maps 

{(No M) ft : p is sufficiently small} 
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is a family of rank one maps to which the theory of [30] and [34] directly applies. In this 
parameter range, the stable and unstable manifolds of the perturbed saddle do not intersect. 
The dynamical properties of the periodically-perturbed equations are determined by the 
magnitude of the forcing frequency u. When the forcing frequency u is small, there exists an 
attracting torus in the extended phase space for all yU sufficiently small. In particular, there 
exists an attracting torus consisting of quasiperiodic solutions for a set of fj, with positive 
Lebesgue density at // = 0. As u; increases, the attracting torus is disintegrates into isolated 
periodic sinks and saddles. Increasing the magnitude of the forcing frequency oj further, the 
phase space is stretched and folded, creating horseshoes and strange attractors. We prove 
in particular that these are strange attractors with SRB measures. SRB measures represent 
visible statistical law in nonuniformly hyperbolic systems. The chaos associated with them 
is both sustained in time and observable. First constructed for uniformly hyperbolic systems 
by Sinai [26], Ruelle [25], and Bowen [7], SRB measures are the measures most compatible 
with volume when the volume is not preserved. See [35] for a review of the theory and 
applications of SRB measures. 

In this paper, we focus exclusively on the scenario of rank one chaos. See Theorems 1 and 
2 in [31] for results concerning the other scenarios described above. We remark that [30] not 
only proves the existence of SRB measures, but also establishes a comprehensive dynamical 
profile for the maps with SRB measures. This profile includes a detailed description of the 
geometric structure of the attractor and statistical properties such as exponential decay of 
correlations. We have opted to limit the statements of our theorems to the existence of SRB 
measures, but all aspects of this larger dynamical profile apply. 

This paper is not only about the generic existence of rank one attractors in periodically- 
forced second order equations. Explicit, verifiable conditions are formulated. Based on the 
theorems of this paper, the first named author has proven the existence of rank one chaos 
in a Duffing equation of the form 

cPq , 2 dq o 

-j^ + ia-bq )—-q + q = fismut 

and in a periodically-forced pendulum of the form 

d 2 6 ^d6 . „ 

— - — o— + sin — a + fi sin ut. 

dt l dt 

These results will be presented in separate papers. 

The analysis in this paper is not sensitive to the particular form we have chosen for the 
forcing. We work with the forcing function Vv,pv because the resulting analysis is relatively 
transparent. A theorem analogous to our main theorem holds for a general class of forcing 
functions. 

This paper is organized as follows. We state our results precisely in Section 2. In Section 3 
we discuss a model of Afraimovich and Shilnikov. Sections 4-7 are devoted to the proof of 
the main theorem. 

D. Acknowledgment. Our method is motivated by a paper of V.S. Afraimovich and L.P. 
Shil'nikov published almost thirty years ago [1]. Afraimovich and Shil'nikov observed that 
for periodically-forced systems with dissipative homoclinic loops, the dissipation around the 
fixed point could potentially put the flow-induced return maps into the category (in our 
terminology) of rank one maps. In this paper we basically start from where they stopped, 
turning an insightful observation into a theorem one can use to analyze concrete equations. 
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We are deeply indebted to Afraimovich for bringing his previous work with Shilnikov [1] to 
our attention. See also [2]. We also thank Kening Lu and Lai-Sang Young for motivating 
conversations related to this work, and particularly Lai-Sang Young for connecting us to 
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2. Statement of Results 

Let (x, y) G M 2 be the phase variables and t be the time. We start with an autonomous 
system 



(2.1) 



^ = -ax + f(x,y) 
^ = Py + g(x,y) 



where / and g are real analytic at (x, y) = (0, 0) and /(0, 0) = g(0, 0) = d x f(0, 0) = 
d y f(0, 0) = d x g(0, 0) = d y g(0, 0) = 0. We assume that a and f3 satisfy a certain Diophantine 
nonresonance condition and that (x,y) = (0,0) is a dissipative saddle point. Namely, we 
assume the following. 

(HI) Nonresonant dissipative saddle. 

(a) There exist di, d<i > such that for all m, n G Z + , we have 

\ma - nj3\ > di(\m\ + \n\y d ' z . 

(b) < (3 < a. 

We also assume that the positive x-side of the local stable manifold of (0, 0) and the positive 
?/-side of the local unstable manifold of (0, 0) are included as part of a homoclinic solution 
which we denote as x = a(t), y = b(t). Let 

£ = {£(t) = (a(t),b(t)) G K 2 : t G R}. 

We further assume that f(x, y) and g(x, y) are C 4 in a sufficiently small neighborhood of I. 

To the right side of equation (2.1) we add a time-periodic term to form a non-autonomous 
system 



(2.2) 



dx 

— = -ax + f(x, y) - /i(ph(x, y) + smut) 
dy 

= Py + g(x, y) + n{ph(x, y) + sin^t) 



where //, p, and u are parameters. We assume that h(x,y) is analytic at (x, y) = (0,0) and 
C 4 in a small neighborhood of the homoclinic loop i. The parameter satisfies ^ C 1 
and controls the magnitude of the forcing term. The prefactor p and the forcing frequency uo 
are much larger parameters, the ranges of which we will make explicit momentarily. Observe 
that the same forcing function is added to the equation for y but subtracted from the equation 
for x. We do this to facilitate the application of our theorem to a certain concrete second 
order system. The analysis in this work is by no means limited to these particular forcing 
functions. 
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To study (2.2), we introduce an angular variable 9 e S 1 and write it as 



( dx 



(2.3) 



We denote 



_ = -ax + f(x, y) - p(ph(x, y) + sin 9) 



dy 
dt 

de 
It 



(3y + g(x, y) + fi(ph(x, y) + sin 9) 



= cu. 



(u(t),v(t)) 



where £(t) = (a(t),b(t)) is the homoclinic loop of equation (2.1). The vector (u(t),v(t)) is a 
unit vector tangent to i at £(t). Define 

E(t) = v 2 (t)(-a + d x f(a(t),b(t))) + u 2 (t)(P + d y g(a(t), b(t))) 

- u(t)v(t)(d y f(a(t),b(t)) + d x g{a{t), b(t))). 



(2.4) 



The quantity E(t) measures the rate of expansion of the solutions of equation (2.1) in the 
direction normal to £ at £{t) (see Section 4.2). In matrix form, we have 



E(t) = {v(t) -u(t)) 



-a + d x f(£(t)) d y f{£{t)) 
d x g(£(t)) (3 + d y g(e(t)) 



v(t) 
-u(t) 



Define 



(2.5) 



A 
C 

s 



(u(s) + v(s))h(a(s),b(s))e-ti E ^ dT ds 
(u(s) + v(s)) cos(oos)e- % e{t) dr ds 
(u(s)+v(s))sm(cus)e-f° EiT)dT ds. 



The integrals A, C, and S are all absolutely convergent (see Lemma 4.3). They describe 
the relative positions of the stable and unstable manifolds of the perturbed saddle. See 
Fig. 4. The quantity pA/j, measures the average distance between the stable and unstable 
manifolds and p{C 2 + S 2 ) 2 measures the magnitude of the oscillation of the unstable manifold 
relative to the stable manifold. 




Fig. 4. The geometric meaning of the integrals A, C, and S. 
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We assume that A, C, and S satisfy the following nondegeneracy conditions. 

(H2) Nondegeneracy conditions on A, C, and S. 

(a) Aj^O. 

(b) C 2 + S 2 ^ 0. 

Given equation (2.2) satisfying (HI) and (H2), we let 



202 VC 2 + S 2 396 VC 2 + S 2 

Pl ~""99" A ' P2 ~~101 A ' 

We also let 

I = {z G R, \z\ < Ku} 
for some K > 1 sufficiently large independent of \x and 

E = {£(0) + (u(0), -u(0))z G M 2 : z G /} x S 1 . 
The following is the main theorem of this paper. 

Theorem 2.1. Assume that (2.2) satisfies (HI) and (H2)(a). There exists ujq > such 
that if uj G R satisfies (H2)(b) and \u\ > uj , then for every p G [^1,^2] we /ia^e i/ie 
following. 

(1) For fi sufficiently small, equation (2.3) induces a well-defined return map 9^ : E — > E. 
There exists a set A WiP of values of \i with positive lower Lebesgue density at fi = 
s^c/i t/iat for every \i G A W((0; 9^ admits a strange attractor that supports an ergodic 
SRB measure v. Furthermore, Lebesgue almost every point on E zs generic with 
respect to v. 




Fig. 5. The Poincare section E. 

We recall that an ^-invariant Borel probability measure v on E is an SRB measure if 
has a positive Lyapunov exponent ^-almost everywhere and if the conditional measures of v 
on unstable manifolds are absolutely continuous with respect to the Riemannian measures 
on these unstable leaves. SRB measures represent visible statistical law in chaotic systems. 

Remark 2.1. As an important condition to be verified, (H2) does not cast doubt on the 
abundance of the type of strange attractor proved to exist in this paper. By properly 
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adjusting the sign of h(x, y) according to the sign of u(s) + v(s) on £, we can easily achieve 
A^O. Hypothesis (H2)(b) requires that the Fourier spectrum of the function 

R( s ) = ( u (s)+v(s))e-f° E ^ dT 

is not identically zero on the frequency range higher than u> . Since R(s) decays exponentially 
as a function of s, the Fourier transform R(£) is analytic in a strip containing the real £-axis 
by the Paley- Wiener theorem. It follows that R(g) = for at most a discrete set of values 
of £ unless R(s) is identically zero. 

3. A MODEL OF AFRAIMOVICH AND SHILNIKOV 

In this section we study a model introduced by Afraimovich and Shilnikov in [1]. See 
also [2]. This simple model allows us to illustrate the steps of the proof of the main theorem 
without needing to deal with technical complexity. The return maps of Afraimovich and 
Shilnikov are derived in Section 3.1. In Section 3.2 we prove that these return maps are rank 
one maps in the sense of [30] and [34]. 

3.1. Derivation of return maps. We begin by describing an unperturbed system of dif- 
ferential equations. Let / : R 2 — > R and g : R 2 — > R be C°° functions and let a, (3 G K 
satisfy < (3 < a. Define 



(3.1) 



^ = -ax + f(x,y) 

% = Py + g(x,y). 

at 



We assume that the functions / and g satisfy f(x,y) = g(x,y) = for all (x,y) G B(0,2e) 
where < e < 1. This means that equation (3.1) is linear in a neighborhood of 0. We also 
assume that equation (3.1) admits a homoclinic solution i = {£(t) : t G R} containing the 
segments {(0, y) : < y < 2e} and {(x, 0) : < x < 2e}. 

Let S 1 = [0, 27r) denote the unit circle and let p, q : R 2 x S 1 — > R be C°° functions such 
that p = q = on B(0, 2e) x S 1 . We now introduce the perturbed system 



(3.2) 



^ = ax + f{x, y) + np(x, y, 9) 
dy 

= Py + g(x, y) + nq(x, y, 9) 

d9 

— = UJ. 

dt 



Here u G 1 is the frequency of the forcing functions and fi > represents the strength of 
the perturbation. We assume that /x and e satisfy ^ /i <ti e < 1. 

The orbit 7 = {(0, 0, 9) : 9 e S 1 } is a hyperbolic periodic orbit of equation (3.2) for all /i. 
For fj, = 0, T = £ x S 1 is the stable manifold and the unstable manifold of 7. We define the 
Poincare sections 

£- = {( x , y, 9) : ^ x ^ dfi, y = e, 9 G S 1 } 

S+ = {(x,y,9):x = e, C 2 ~V < V < C 2f i, 9 G S 1 } 

where \x G [0,/io], C\ > is such that Ci/io e, and C2 is suitably chosen. We study 
a situation in which one can define flow-induced maps M : £~ — > S + and : S + — > E _ 
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(see Section IB and Fig. 3). The composition HoM produces a one-parameter family 
{3^ = N o M : /i G [0, /i ]} of maps from £~ to £~. 

The map 3V : £+ — ► S _ . The flow from S + to £~ is defined by the differential equations 

(3.3a) — = —ax 

v 1 dt 

(3.3b) *=ft, 

,,,, <*> 
(3.3c) — = a;. 

Let (e,y,9) G S + . Let denote the time at which the orbit emanating from (e,y,6) 
intersects £~. Write N(e,y,Q) = (xi,e,9i). Integrating (3.3b), we have e = e^ T ^y, so 



a a 



T(y) = 4log(£y r ). Integrating (3.3a) yields x\ — e aT ^e = e 1 PyP . The local map N is 



P 

therefore given by 



(3.4) 



Xi — e Pyl 3 

e^e + ^iogiey- 1 )- 



The map M : £ — > £+. Let (x ,e,9 ) G X . Write M(x ,£, #o) = (z,y,9)- We assume 
that for /i G [0, //o], 

= 0o + fi + /#(z o ,0o)- 

Here < A < 1 and £i > are fixed. The functions </? and ip are C°° functions on £~. We 
assume that <p(xo, 9q) > for all (xo, 9q) G This ensures that the second scenario of Fig. 
1, namely the scenario in which the stable and unstable manifolds are pulled apart by the 
periodic forcing. More precisely, we assume that p and q are such that 

ip(x ,9 ) = £ 2 

ip(x ,9 ) = B(l + Asm9 ). 
Here £ 2 G K, B > 0, and < A < 1. The global map M is therefore given by 

y = \x + fiB(l + A sin 9 ) 
9 = 9 + £i + /x&. 



(3.5) 



Let us not worry about the viability of these assumptions, enduring for the moment the 
possibility that, at worst, no differential equation satisfies all of our assumptions. 

The map J (1 = NoM: S _ — > S _ . Let (x ,e,9 ) G £~. Computing ^^(x 0: e,9o) = 
(xi,e,9i) using (3.4) and (3.5), we have 

Xl = [Ax + iiB{\ + Asin(# )] f 
0i = #o + 6 + A^ 2 + % log 



/3 VAzo + //£(l + Asin(0 o )) 
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Using the spatial rescaling x i— > fiX, we obtain 

CX.CH. CX 

(3.6a) X 1 = £ 1 ~Pfj,P~ 1 [\X + £(l + Asin(0 o )F 

(3.6b) X = 0o + £i + + £ log f — ^ . ) . 

fj \XX + B{1 + Asm(6 )) J 

Using this formula for 3^, Afraimovich and Shil'nikov [1, 2] conclude that 3^ has a horse- 
shoe for large a;. 

3.2. Theory of rank one attractors. In this subsection we first introduce admissible rank 
one maps following [34] and we then prove that {3^} is an admissible family of rank one 
maps using the techniques of [31]. 

A. Misuirewicz maps and admissible ID families. The definition of an admissible family 
of ID maps is rather long and technical. It could therefore present a nontrivial hurdle for the 
reader. We feel obligated to present this definition for completeness. Readers wishing to skip 
the material on admissible ID families can safely jump to Proposition 3.1. Proposition 3.1 
contains the only result from the ID aspect of rank one theory that we need for the results 
of this paper. 

We start with Misuirewicz maps. For / G C7 2 (S 1 ,S 1 ), let C = C(f) = {/' = 0} denote 
the critical set of / and let Cs denote the (^-neighborhood of C in S 1 . For x G S 1 , let 
d(x, C) = minfgc \x — x\. 

Definition 3.1. We say that / G C^S^S 1 ) is a Misuirewicz map and we write / G £ if the 
following hold for some Sq > 0. 

(1) Outside of Cs - There exist Ao > 0, M G Z + , and < c < 1 such that 

(a) for all n > M , if x, f(x), ■■■ , f n ~\x) C So , then |(/ n )'(a0| > e A °"; 

(b) if x, f(x), • • • , f n ~\x) C So and f n (x) G C So for any n, then |(/")'(^)| > c e A «". 

(2) Inside Cs - 

(a) We have f"{x) ^ for all x G Cg . 

(b) For all x G C and n > 0, d(f n (x), C) > S . 

(c) For all x G C$ \ C, there exists po(x) > such that f j (x) G" Cs for all j < po(x) 
and |(/ Po(:c) )'(a;)| > C q 'ei^W. 

We remark that Misurewicz maps are among the simplest maps with nonuniform ex- 
pansion. The phase space is divided into two regions, Cs and S 1 \ Cs Q - Condition (1) in 
Definition 3.1 says that on 8 1 \Cg , f is essentially uniformly expanding. Condition (2c) says 
that for x G C$ \ C, even though is small, the orbit of x does not return to Cg again 

until its derivative has regained a definite amount of exponential growth. In particular, if n 
is the first return time of x G Cs to Cs , then \{f n )'{x)\ > CQ 1 e~3 Xon . 

We now define admissible families of ID maps. Let F : S 1 x [a 1: a 2 ] -> S 1 be a 
C 2 map. The map F defines a one-parameter family {f a G C 2 (§ 1 ,§ 1 ) : a G [ai,a 2 ]} via 
fa( x ) — F(x,a). We assume that there exists a* G (ai,a 2 ) such that f a * G £. For each 
c G C(/ a »), there exists a continuation c(a) G C(/ a ) provided a is sufficiently close to a*. 

Let C(/ a *) = {^'(a*),...,^^)}, where c W(o*) < c( m )(a*) for 1 ^ % ^ g - 1. For 
c(a*) G C(f a *), we define /5(o*) = / a *(c(a*)). For all parameters a sufficiently close to a*, 
there exists a unique continuation /3(a) of /3(a*) such that the orbits 

{fS*((3(a*)) : n ^ 0} and {/^(a)) : n > 0} 
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have the same itineraries with respect to the partitions of S 1 induced by C(/ *) and C(f a ). 
This means that for all n > 0, j™(0(a*)) G (c^(a*),c^ +1 \a*)) if and only if f a l (P(a)) G 
(c^(o), c^ +1 ^(a)) (here c*- 9 " 1 " 1 - 1 = c^). Moreover, the map a i— > /3(a) is differentiable (see 
Proposition 4.1 in [33]). 

Definition 3.2. Let F : S 1 x [ai,a 2 ] — > S 1 be a C 2 map. The associated one-parameter 
family {f a : a G [01,02]} is admissible if 

(1) there exists a* G (cti, a 2 ) such that f a * G £; 

(2) for all c G C(/ *), we have 



(3-7) £(c)=^(/ a (c(a))-/?(a)) 



7^0. 



The next proposition contains all that we need from the ID aspect of rank one theory for 
this paper. 

Proposition 3.1 ([31, 32, 17]). Let : S 1 — > R fre a C 3 function with non- degenerate 
critical points and let $(#, a) : S 1 x [a , ai] — > R 6e swcn £nat 

ll^(^,a)l|c3(Six[ao,oi]) < 



100 

We define a one parameter family of circle maps {f a : a G [0, 27r]} fry 

/ o (0) = + 0(0,a) + a + K*(0) 

where X is a constant. There exists K , determined by \& alone, such that if % > K , then 
{fa} is an admissible family of ID maps. 

The special case of this proposition in which $(0, a) = was first proved in [31]. That 
proof can easily be extended to prove Proposition 3.1. See also Proposition 2.1 in [32] and 
Appendix C in [17]. 

B. Admissible families of rank one maps. We now move to the 2D part of the setting 
of [30] and [34]. Let / be an interval. Let B Q C R be a set with a limit point at 0. A 
2-parameter C 3 family {F a t b(X, 9) : a G [ao,ai], b G -Bo} of 2D diffoemorphisms defined on 
S = / x S 1 is an admissible rank one family if the following hold. 

(CI) There exists a C 2 function F a> o(X, 6) of (a, X, 6) such that, as b — > 0, 

\\F a ,b(X,9) - (0, F afi (X, #))||c 3 ([a ,ai]x£) ~ > 0. 

(C2) {f a {9) = F a0 (0, 9) : a G [a , a±]} is an admissible ID family. 

(C3) For all a G [a , 01], at the critical points of the ID map f a (9) we have 



d 

-r^F a ,o(X,9) 



7^0. 

x=o 



The following is the main result of [30] and [34] for a given admissible rank one family F a ^ 
of 2D maps. 

Proposition 3.2 ([30, 34]). Let F a> b : X — > S fre an admissible rank one family. There exists 
b > saca i/iai /or a// |6| < 6, inere exists a set 0/ values of a with positive Lebesgue 
measure such that for a E Ab, F a ^ admits an ergodic SRB measure v. If we also have 
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A > In 10, where A is as in Definition 3.1, then v is the only ergodic SRB measure 1 that 
F ajb admits on E. 

More is true if the global distortion bound (C4) holds. 

(C4) There exists C > such that for all a G [a , ai], b G B , (X, 9) G E, and (X', 9') G E, 
we have 

detDF atb (X,6) 



det DF a>b (X>,9>) 



< C. 



Proposition 3.3 ([30]). Let F a ^ be an admissible rank one family satisfying (C4) and 
suppose that X > In 10, where Ao is as in Definition 3.1. Then for all \b\ < b and a G A b , 
Lebesgue almost every point in E is generic with respect to the unique ergodic SRB measure 
on E. 

C. {3^} is an admissible family of rank one maps. We show that {3^} satisfies 
the hypotheses (C1)-(C4). Letting fx — > in (3.6a), we see that X x — > because a > (3. 
However, the term % \og{fi~ 1 ) (mod 2ir) fails to converge as fx — > 0. The fact that 9\ is 
computed modulo 2n allows us to introduce the parameter a and thereby obtain a two- 
parameter family {F ab } with a well-defined ID singular limit. 

We regard p = log(/x _1 ) as the fundamental parameter associated with {3^}. Notice that 
we now have p G [log(/XQ x ), 00). Think of // = e~ p as a function of p. Define 7 : (0, fx ] — > K. 
by 

Let iVeN satisfy |log(^g 1 ) < iV. Let (/i n ) be the decreasing sequence of values of fx such 
that 7(/i„) = N + 2n{n — 1) for every n G N. We think of fx as a measure of dissipation and 
we therefore set b n = fx n . For a G S 1 and n£N, define 

fx(n,a) = 7 _1 (7(^n) + a) 

p(n, a) = log(/i(n, a) -1 ) = log^" 1 ) + -a. 

The map F afin is defined by F aK = 3p( n , a ). 

The family {-F a ,fc„} has a well-defined singular limit. As n — > 00, F a bn converges in the C 3 
topology to the map F a defined by 

f$(x o ,0o) = o 

F a 5(Xo, 6 ) = 9 + 6 + ^ log(e) + a - ^ log(AX + S(l + Asin(0 o ))). 
This proves (CI). 

Restricting Fj 2 J to the circle {(X ,9 ) : X = 0}, we obtain the one-parameter family of 
circle maps 

f a {9) = 9 + ^1 + ^ log(e) + a - ^ log(S(l + A sin(0))). 

It follows directly from Proposition 3.1 that f a is an admissible family of ID maps provided 
is sufficiently large. This proves (C2). Hypotheses (C3) and (C4) follow from direct 
computation. 



1 This is proved in [31]. 
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We have shown that the family {F abn } is an admissible rank one family and therefore 
Propositions 3.2 and 3.3 apply. We conclude that if \u>\ is sufficiently large, then there exists 
a set of positive Lebesgue measure such that for p e A w , admits a strange attractor 
on £~ with an ergodic SRB measure v and Lebesgue almost every point on X~ is generic 
with respect to v. Furthermore, the set A w has positive lower Lebesgue density at 0, meaning 
that 

where | • | denotes Lebesgue measure. 

4. Standard forms around the homoclinic loop 

In this section we introduce a sequence of coordinate changes to transform equation (2.3) 
into certain standard forms. In Section 4.1 we work in a sufficiently small neighborhood 
U £ of (0,0) in the (x, y) -plane. In Section 4.2 we work in a small neighborhood around 
the entire length of the homoclinic loop £ outside of U-yi. In Section 4.3 we define the 
Poincare sections which we will use to compute the flow-induced maps. Points on S 1 * 1 are 
represented differently by various sets of variables introduced in Sections 4.1 and 4.2. We 
discuss the issue of coordinate conversion in Section 4.3. 

In the rest of this paper, a, /3, p <E [pi, P2} and co > oo (it suffices to assume u > 0) are 
all regarded as fixed constants. The size of the neighborhood on which all of the coordinate 
transformations in Section 4.1 are performed is determined by a small number e > 0. The 
quantity e is also regarded as a fixed constant. We regard p as the only parameter of 
equation (2.3). 

Two small scales. The quantities |i C £ < 1 represent two small scales of different 
magnitude. The quantity e represents the size of a small neighborhood of (x, y) = (0, 0) in 
which the local analysis of Section 4.1 is valid. Define 

U £ = {(x, y):x 2 + y 2 < 4e 2 } and U £ = U e x S 1 . 

Let L + and — L~ be the respective times at which the homoclinic solution £{t) enters Ui £ in 
the positive and negative directions. The quantities L + and L~ are completely determined by 
e and £. The parameter p (p <C e) controls the magnitude of the time-periodic perturbation. 

Notation. Quantities that are independent of phase variables, time and p are regarded as 
constants and K is used to denote a generic constant, the precise value of which is allowed 
to change from line to line. On occasion, a specific constant is used in different places. We 
use subscripts to denote such constants as K , K 1} - ■ ■ . We will also distinguish between 
constants that depend on e and those that do not by making such dependencies explicit. A 
constant that depends on e is written as K{e). A constant written as K is independent of e. 

4.1. Standard form near the fixed point. In this subsection we study equation (2.3) 
in a sufficiently small neighborhood of (0, 0) in the (x, y)-plane. We introduce a sequence 
of coordinate changes to transform equation (2.3) into a certain standard form. Table 1 
summarizes the purpose of each coordinate transformation. 



A. First coordinate change: (x,y) — > (£,77). Let (£,r/) be such that 
(4.1) € = x + qi(x,v), V = V + Q2(x,y) 
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Table 1. Transformations near the fixed point. 



Tr ansf or mat ion 


Purpose 


(x,y) -> (£,ri) 


linearize the flow defined by (2.1) in a neighborhood of (0,0) 


(£, V )^(X,Y) 


standardize the location of the hyperbolic periodic orbit 


(X,Y)^(X,Y) 


flatten the local invariant manifolds 


(X,Y)-(X,Y) 


rescale by the factor p~ l 



where qi(x, y) and q2(x, y) are analytic terms of order at least two in x and y. Formula (4.1) 
defines a near-identity coordinate transformation (x,y) — > (£,?/), the inverse of which we 
write as 

(4.2) x = ^ + Q 1 ^,r t ), y = r ] + Q 2 (Z, V ). 

Proposition 4.1. Assume that a and (3 satisfy the nonresonance condition (HI )(a). Then 
there exists a neighborhood U of (0,0), the size of which is completely determined by equa- 
tion (2.1) and d\ and d 2 in (Hl)(a), such that on U there exists an analytic coordinate 
transformation (4-1) that transforms equation (2.1) into the linear system 

d£ _ dq 

dt ~ dt 

Proof. See [12] for a proof. ■ 

We now use the coordinate transformation of Proposition 4.1 to transform equation (2.3). 
Observe that by definition, qi(x,y) and q 2 (x,y) satisfy 

(4.3a) (1 + d x qi(x, y))(-ax + f(x, y)) + dyq^x, y)(f3y + g(x, y)) = -a£ 

(4.3b) (1 + d y q 2 (x, y))((3y + g(x, y)) + d x q 2 (x, y)(-ax + f(x, y)) = fir). 

We derive the form of (2.3) in terms of £ and rj. We have 

^ = (\ + d x q 1 (x,y))(-ax + f(x,y) - p,(ph(x,y) + sin#)) 

+ d y qi{x,y)(Py + g(x,y)+ii(ph(x,y) + sm6)) 
= -a£ - n(l + d x qi(x,y) - d y q 1 (x,y))(ph{x,y) + sintf) 
where (4.3a) is used for the second equality. Similarly, we have 
dri 

— = (1 + d y q 2 (x,y))((3y + g(x,y) + p(ph(x,y) + sin 6)) 

+ d x q 2 (x, y)(-ax + f(x, y) - p(ph(x,y) + sin#)) 

= Pv + M 1 + d yl2(x,y) - d x q 2 (x,y))(ph(x,y) + sin6>). 

Writing the functions of x and y as functions of £ and n using (4.2), the form of (2.3) in 
terms of £ and rj is given by 

"37 = -«? 



(4.4) 



dt 
drj 

lit 

de 

~dt 



/i(l + M^))G*#(^) + sin0) 
P V + p(l + h 2 ^,rj))(pH(C,v) + sine) 



= uo 
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where hxfarj) = d x qi{x,y) - d v qi(x,y), h 2 (£,rf) = d y q 2 (x,y) - d x q 2 (x,y) are such that 
^(0,0) = h 2 (0,0) = and H(^,rf) = h(x,y). 

B. Second coordinate change: (£,77) — > (X,Y). With the forcing added, the hyperbolic 
fixed point (x,y) = (0,0) of equation (2.1) is perturbed to become a hyperbolic periodic 
solution of (2.2) period 2tiuj~ 1 . We denote this periodic solution in (£, 77, 6>)-coordinates as 
£ = p4>(9; p), 77 = ptp(9;p). 

Proposition 4.2. For equation (4-4), there exists a unique solution of the form 

£ = p<j)(9;p), 77 = /i^(0;/i), = ^rj 

satisfying 

<j>(6; p) = <f)(9 + 27r; /i), t/;(0; A*) = ^(0 + 2tt; /x). 

TTie C 3 norms of the functions <f>(6; p) and ifi(9;p), regarded as functions of 9 and p, are 
bounded by a constant K . 

Proof. Write = 0(0; p), ip — ip(9; p). The functions and -0 should satisfy 

< 4 - 5 > 2 

U l§ = P^ + ( 1 + MM v4>))(pH(p4>, + sin 0). 
From (4.5) it follows that 

0(0; p) = e-™- 1 ^^ fi) - cu' 1 f e™- 1 ^ [1 + h^s; p),p^{s; fi))] x 

[pH(fi(j){s] fi), pip(s; fi)) + sin s] ds 
iP(9; fi) = e^'^V^o; fi) + co' 1 f e-^ s ~ e \\ + h 2 (p(f>(s; p),p^{s; fi))] x 

JOo 

[pH(fi(j)(s; fi), fiifj(s; fi)) + sin s] ds. 

To solve for and ip we let 9 = O + 2n and set 0(0o + 27r; /j) = 0(0o; fi), ip(9o + 2n; fi) = 
ip(9 ; fi), obtaining 

*(M = i / e a ^ s - 2 ^[l + h 1 (fi4>(s + 9;fi),fii;( S + 9;fi))]x 

1 — e Jo 

[pH(fuf)(s + 9; fi),fii{j(s + 9; fi)) + sin(s + 9)] ds 

( 4 - 6 ) 1 ,2. 

W>ri = T=&£F* J o e-^ s - 2 ^[l + h 2 (fi ( f ) (s + 9 ] fi),fi^ S + 9 ] fi))]x 

[pH{fi(f){s + 0; p),pip(s + 9; p)) + sin(s + 0)] ds. 

The existence and uniqueness of 0(0; p) and ^(0; p) follows directly from an application of 
the contraction mapping theorem to (4.6). The asserted bound on partial derivatives with 
respect to 9 and p follows from differentiating (4.6) with respect to 9 and p. ■ 

We now introduce new variables (X, Y) by defining 
(4.7) X = £-p(f>(9;p), Y = V -p^(9;p). 
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We have 

= - a X - ap<p - fiuj^r - p(l + h(X + pep, Y + pip))(pH(X + pcp,Y + flip) + sin#) 
at ciu 

dY dip 

— = pY + (5pip - fitu-j- + + h 2 (X + fjuf>, Y + pip))(pH(X + p<p,Y + pip) + sin0). 
dt du 

Using (4.5), the form of (2.3) in terms of X, Y and 9 is given by 

( dX 



(4.8) 



df aX I pF(X.Y9:p) 

— — to 
dt 



where 

F(X, y, 9; p) = -[hiX + p<P,Y + flip) - h^, piP)](pH(X + p<f>, Y + [lib) + sin 9) 
- p(l + /ii(/i0, pip))(H(X + pcp,Y + pip) - H(fjuf>, pip)) 

G(X,Y,9;p) = [h 2 (X + pcp,Y + pip) — h 2 (pcp, pip)](pH(X + pep, Y + pip) + sin 9) 
+ p(l + h 2 (p(P, pip))(H(X + p<P,Y + pip) - H(p<P, pi>)) 

are such that F(0, 0, 9; p) = G(0, 0, 9; p) = 0. Observe that in the new coordinates (X, Y, 9), 
the solution £ = p<p(9; p), rj = pip {6; p) is represented by X = Y = 0. We remark that on 

{(X,Y,9;p):\\(X,Y)\\<e, 9eS\ 0^p^p }, 

(1) F(X, Y, 9; p) and G(X, Y, 9; p) are analytic functions bounded by Ke; 

(2) it follows from Proposition 4.2 that the C 3 norms of both F and G as functions of 
(X, Y, 9) and p are bounded by a constant X. 

C. Third coordinate change: (X, Y) -> (X, Y). The periodic solution (X, Y, 0) = (0, 0, cut) 
of equation (4.8) has a local unstable manifold, which we write as 

X = pW u (Y,9;p), 

and a local stable manifold, which we write as 

Y = pW s (X,9;p). 

Proposition 4.3. There exists e > and p = po(e) > such that W U (Y,9; p) and 
W S (X, 9; p) are analytically defined on 

(-£,£) X S 1 X [0,p ] 

and satisfy 

W u (0, 9; p) = 0, W s (0,9;p) = 0. 

TTie C 3 norms of W U (Y,9; p) and W s (X,9;p), regarded as functions of all three of their 
arguments, are bounded by a constant K. 

Proof. We regard X, Y, 9, and p in equation (4.8) as complex variables. The existence and 
smoothness of local stable and unstable manifolds follows from the standard argument based 
on the contraction mapping theorem. See [12] for instance. ■ 
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By definition, W U (Y, 9; fi) satisfies 

-aW u (Y, 9; fi) + F(fiW u (Y, 9; fi),Y, 9; fi) = ud e W u (Y, 9; fi) 
+ d Y W u (Y, 9; n){0Y + fiG(fiW u (Y, 9; fi), Y, 9; fi)). 
Similarly, W S (X, 9; fi) satisfies 

f3W s (X, 9; fi) + G(X, fiW s (X, 9; fi), 9; fi) = ud e W s (X, 9; fi) 
+ d x W s (X, 9; fj)(-aX + fiF(X, fiW s {X, 9; fi),9; fi)). 
Define the new variables X and Y by 
(4.11) X = X - fiW u (Y,9;fi), Y = Y — fiW s (X, 9; fi). 

By using (4.8), (4.9), and (4.10), the form of (2.3) in terms of (X, Y,9) is given by 



(4.9) 



(4.10) 



( ^ = (-a + fiF(X,Y,9;fi))X 



(4.12) 



dY 

~dt 
d9 

~db 



(/3 + /iG(X,Y,0;/i))Y 



cu. 



where F and G are analytic functions of X, Y, 9, and fi defined on U e x S 1 x [0,/xo]- The 
C 3 norms of F and G are bounded by a constant K. Tracing back to the variables (£,77), 
we have 



(4.13a) 
(4.13b) 



X = £ - fi (0(0; fi) + W u (rj - rt(6; fi), 9; fi)) 
Y = 7] — fj, (m fi) + W S (C- nm fi), 9; fi)) . 



D. Fourth coordinate change: (X, Y) — > (X, Y). The final coordinate change is a rescaling 
of X and Y by the factor fi -1 . Let 

(4.14) X = /i _1 X, Y = fT 1 Y. 

We write equation (4.12) in X and Y as 

d%. 

= (-a + //F(X,Y,0;//))X 



(4.15) 



dt 
d¥ 

— = (P + fiG(X,Y,9;iu))Y 
d9 



dt 



= uo 



where 

F(X, Y, 9; fi) = F(/iX, fiY, 9; fi), G(X, Y, 9; fi) = G(fjX, fiY, 9; fi) 
are analytic functions of X, Y, 9, and fi defined on 

D = {(X, Y,9, fi) : fi e [0, fi ], (X, Y, 9) e UJ 

where 

U £ = {(X, Y, 9) : || (X, Y) || < 2efi-\ 9 e S 1 }. 

Remark 4.1. We remind the reader that all constants represented by K in Section 4.1 are 
independent of e and fi. 
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4.2. A standard form around the homoclinic loop. In this subsection we derive a 
standard form for equation (2.3) around the homoclinic loop of equation (2.1) outside of 
Wi e 2. Some elementary estimates are also included. 

A. Derivation of equations. Let us regard t in £(t) = (a(t),b(t)) not as time, but as a 
parameter that parametrizes the curve £ in (x, y)-space. We replace t by s and write this 
homoclinic loop as £(s) = (a(s), b(s)). We have 

da(s) 



ds 
db(s) 



ds 

(u{s),v(s)) 



= -aa(s) + f(a(s),b(s)) 
= Pb(s)+g(a(s),b(s)). 



d_ 

ds 



£(s) 



- 1 d_ 
ds 



1(8). 



u(s) 
v(s) 



i(s) + f(a(s),b(s)) 



(4.16) 
Define 

We have 

(4.17) 

Let 

e( s ) = (v(s),-u(s)). 

The vector e(s) is the inward unit normal vector to £ at £(s). We now introduce the new 
variable z such that 

(x,y) = £(s) + ze(s). 

That is, 

(4.18) x = x(s, z) = a(s) + v(s)z, y = y(s, z) = b(s) — u(s)z. 

We derive the form of (2.3) in terms of the new variables (s, z) defined through (4.18). 
Differentiating (4.18), we obtain 



y/(-aa(s) + f(a(s), b(s)))* + W + 9(a(s)b(sW ' 

Pb(s)+g(a(s),b(s)) 
^(-aa(s) + f(a(s), b(sW + W*) + 9(a(s), b(s)))* 



(4.19) 



where u'(s) 



^ = (-aa(s) + f(a(s),b(s)) + v'(s)z)^ + v(s)^ 

f f = (Pb(s) + g(a(s),b(s)) ~ u'(s)z)^ - u(s)^ 
and v'(s) = Denote 



F(s, z) = —a(a(s) + zv(s)) + f(a(s) + zv(s), b(s) — zu(s)) 
G(s, z) = f3(b(s) — zu(s)) + g(a(s) + zv(s), b(s) — zu(s)) 
H(s, z) = h(a(s) + zv(s), b(s) — zu(s)). 
Using (2.3) and (4.19), we have 

ds _ v(s)G(s, z) + u(s)F(s, z) + fJ,(v(s) - u(s))(pEL(s, z) + sin#) 

dt ~ 

dz 



y/F(s, 0) 2 + G(s, 0) 2 + z(u(s)v'(s) - v(s)u'(s)) 
= v(s)F(s, z) - u(s)G(s, z) - fi(u(s) + v(s))(pM(s, z) + sin#). 
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We rewrite these equations as 

ds n , n \ V(v(s) -u(s))(pH(s,O) + sin0) 
— = 1 + zw^s, z, 9; /i + , = 

( 42 °) ^ = £(s)z + 2 2 w 2 (s, -2) - /i(u(s) + u(s))(pM(s, z) + sin 0) 



de 
It 



where 



E( s ) = v 2 {s){-a + d x f(a(s), b(s))) + u 2 (s)((3 + d y g(a(s),b(s))) 

- u(s)v(s)(d y f(a(s),b(s)) + d x g(a(s),b(s))) 

H(s,0) = h(a(s),b(s)). 

Equation (4.20) is defined on 

{s G [-2L",2L + ], fie [0,H, G S 1 , |z| < K {e)fi}, 

where K (e) is independent of ji. The C 3 norms of the functions Wi(s, z, 0; /i) and w 2 (s, z) 
are bounded by a constant K(e). 

Finally, we rescale the variable z by letting 

(4.21) Z = [i~ x z. 

We arrive at the equations 

ds 

(4.22a) — = l + f iw 1 (s,Z,9; f i) 

dZ 

(4.22b) — = E(s)Z + fiw 2 (s,Z,6;n) - (u(s) + v(s))(pE.(s, 0) + sin0) 

d0 

(4.22c) - = u 

defined on 

T> = {(s,Z,9;ti,):se[-2L-,2L + ], \Z\ ^ K (e), 9 G S\ /x G [0,// ]}. 
We assume that yU is sufficiently small so that 

/.« min (F( S ,0) 2 + G( S ,0) 2 ). 

se[-2L-,2L+] 

The C 3 norms of the functions w\ and W2 are bounded by a constant AT(e) on D. 

System (4.22a)-(4.22c) is the one we need. The function E(s) appears in the integrals A, 
C, and S in (H2). 

Remark 4.2. Observe that all of the generic constants that have appeared thus far in this 
subsection have the form K(e). 

B. Technical estimates. We adopt the following conventions in comparing the magnitude 
of two functions f(s) and g(s). We write f(s) -< g(s) if there exists K > independent of s 
such that \f(s)\ < K\g(s)\ as s — > oo (or — oo). We write f(s) ~ g(s) if in addition we have 
\f(s)\ > K~ 1 \g(s)\. We also write f(s) xs g(s) if 
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as s — > oo (or — oo). 

Recall that £(s) = (a(s),b(s)) is the homoclinic solution for the hyperbolic fixed point 
(0,0) of equation (2.1). The vector (u(s),v(s)) is the unit tangent vector to £ at £(s). 

Lemma 4.1. As s — > +oo ; we have 

(1) a(s) ~ e~ as , a(-s) -< e~ 2fis 

(2) b(s) -< e~ 2as , b(-s) ~ e"^ 
($) u(s) « -1, m(-s) -< e _/3a 

u(s) -< e~ as ; u(-s) « 1. 

Proof. We are simply restating the fact that £(s) — > (0, 0) with an exponential rate —a in 
the positive s-direction along the x-axis and with an exponential rate (3 in the negative 
s-direction along the y-axis. ■ 

Lemma 4.2. Let E(s) be as in (2.4)- As V s1 — > +oo ; we have 

(a) f_V(£(s) + a) ds -< 1 

(b) (E(s) - 0) ds * 1 

(c) J'_ L _E(s)dsm—aL~ 

(d) j Q L+ E{s)ds^ (3L+. 

Proof. Statements (a) and (b) claim that the integrals are convergent as — > oo. For (a), 
we observe that by adding a to E(s), we obtain E(s) + a as a collection of terms, each of 
which decays exponentially as s — > — oo by Lemma 4.1. Similarly, taking /? away from -E'(s), 
we obtain i?(s) — (3 as a collection of terms, each of which decays exponentially as s — > oo. 
For (c) and (d) we write 



/o /-O 
£( s ) ds = -aL~ + / (E(s) + a) ds 
-L- J-L- 
nL+ nL+ 

/ £(s) ds = (3L + + / (E(s)-(3)ds. 
Jo Jo 



Statements (c) and (d) now follow from (a) and (b), respectively. ■ 
Lemma 4.3. All of the integrals defined in (2.5) are absolutely convergent. 
Proof. Let us write 

A= j L ° {u{s) + v{s))h{a{s),b{s))e-^ E{T)dT ds 

J — oo 

+ [ L ° (u(s) + v(s))h(a(s),b(s))e-f° E(T)dT ds 

J-L 

POO 

+ / (u(s) +v(s))h(a(s),b(s))e-^ E{T)dT ds. 

Jl 

We write the first integral as 

/ L \u(s) + v(s))h(a(s),b(s))e as e- f»( E(T ^ dT ds 

J — oo 
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and make L sufficiently large so that \E{t) + a\ < \a for all r G (—00, —L ). This integral 
is convergent since the integrand is < Ke^ as for all s G (—00, —L ). For the convergence of 
the third integral we rewrite it as 

POO 

/ (u(s) +v(s))h(a(s),b(s))e-P s e-fo ( > EiT) ^ )dT ds 

J La 



'L 

and observe that \E(t) — (3\ < f for r G [L , 00) provided that L is sufficiently large. The 
proofs for C and S are similar. 



4.3. Poincare sections and conversion of coordinates. In this subsection we introduce 
the Poincare sections H ± . Since various sets of phase variables have appeared in Sections 4.1 
and 4.2, we also need to know how to explicitly convert coordinates from one set to another 
on T, ± . 

A. The Poincare Sections Y.^. Recall that {£(s) : s G (—00,00)} is the homoclinic loop 
of equation (2.1). Given e > sufficiently small, let L + and — L~ be such that 

C(-L-) = a(-L-) + qi (a(-L-),b{-L-)) = 
rj(-L-) = b(-L-) + q 2 {a{-L-)M-L-)) = e 
{ • ' aL + ) = a(L + ) + qi (a(L + ),b(L + ))=e 

V (L + ) = b(L + ) + q 2 (a(L + ), b(L + )) = 
where ^ and 77 are the variables defined through (4.1). Let 

K = max 

ju6[0,/uo] 

where <j>(9; /i) and ip(9; /i) are as in Section 4. IB. We define two sections in U £ , denoted £~ 
and S + , as follows. 

E" = {(x, y,6):s = -L-, \z\ <: {K + l)/i, 9 G S 1 } 
(4.24) Z+ = {(x,y,6):s = L+, ±(- P A)(K + l)e^ L+ f, < z 

^ lO(-pA)(K + l)e 2/3L+ /x, 9 G S 1 } 

where s and z are as in (4.18). We construct the flow- induced map 3^ in two steps. 

(1) Starting from E _ , the solutions of equation (2.3) move out of U £ , following the homoclinic 
loop of equation (2.1) to eventually hit E + . This defines a flow-induced map from E~ 
to E+, which we denote as M : E~ -> E+. We will prove that M(E~) C E+. 

(2) Starting from E + , the solutions of equation (2.3) stay inside of lt £ , carrying E + into E _ . 
This map we denote as N. 

We define ^ = ^0 M. Observe that the variables (s, Z, 9) of Section 4.2 are suitable for 
computing M and (X, Y, 9) are suitable for computing N. To properly compose 3sf and M, 
we need to know how to convert from (s, Z, 9) to (X, Y, 9) on E 1 * 1 and vice-versa. 
The new parameter p. As stated earlier, we regard /j as the only parameter of system (2.3). 
We make a coordinate change on this parameter by letting p = In jj, and we regard p, not 
/i, as our bottom-line parameter. In other words, we regard p as a shorthand for e p and 
all functions of p are thought of as functions of p. Observe that p G (0, po] corresponds 



22 



QIUDONG WANG AND WILLIAM OTT 



to p G (— oo,ln/i ]. This is a very important conceptual point because by regarding a 
function F(a) of [i as a function of p, we have 

a p F( f i)= l ti ll F(ji). 

Therefore, thinking of F(a) as a function of p produces a C 3 norm that is completely different 
from the one obtained by thinking of F(u) as a function of a. 

Notation 4.1. In order to apply the theory of rank one maps [30, 34], we need to control 
the C 3 norm of 9^. In particular, we must estimate the C 3 norms of certain quantities with 
respect to various sets of variables on relevant domains. The derivation of the flow-induced 
maps {3^} involves a composition of maps and multiple coordinate changes. To facilitate 
the presentation, from this point on we adopt specific conventions for indicating controls 
on magnitude. For a given constant, we write 0(1), O(e), or O(a) to indicate that the 
magnitude of the constant is bounded by K, Ke, or K(e)a, respectively. For a function 
of a set V of variables on a specific domain, we write Oy(l),C?y(e) or Oy{^) to indicate 
that the C 3 norm of the function on the specified domain is bounded by K, Ke, or K(e)a, 
respectively. We choose to specify the domain in the surrounding text rather than explicitly 
involving it in the notation. For example, Ox ,y ,o,h( £ ) represents a function of X , Y , 9, 
and fi, the C 3 norm of which is bounded above by Ke on a domain explicitly given in the 
surrounding text. Similarly, Oz,e, P {u) represents a function of Z, 9, and p, the C 3 norm of 
which is bounded above by K(e)ji. 

B. Conversion on E~. The section S~ is defined by s = —L~. A point q G S~ is uniquely 
determined by a pair (Z, 9). First we compute the coordinates X and Y for a point given in 
(Z, ^-coordinates on S~. Recall that p = hiyu. 

Proposition 4.4. For [i e (0, a ] and (Z,9) G we have 

X = (1 + O e , p (e) + nOz,eM)Z - O e , p {l) 

Y = f i- 1 e + ZAp (l). 

Proof. By definition, s = —L~ on Let q G S~ be represented by (z, 9). Using (4.23), we 
have 

(4.25) ' V ' 

b(-L-) =e + Q 2 (0,e) = e + 0(e 2 ). 

We also have 

(4.26) u(-L-)=0(e), v(-L~) = 1 - 0(e). 
We compute values of X and Y for q. Using (4.23) and (4.25), 

£ = a(-XT) + v(-L~)z + <?i(a(-L - ) + u(-L~)z, b(-L~) - u(-L~)z) 
= v(-L~)z + gi(a(-L") + v(-L~)z, b(-L~) - u{-L~)z) - qi (a(- L ~), K~ L ~)) 
= (l + 0(e) + zhz(z))z. 
Similarly, we have 

rj = b(-L~) - u(-L~)z + q 2 (a(-L~) + v(-L~)z, b(-L~) - u(-L~)z) 
= e - u(-L~)z + q 2 (a(-L-) + v{-L~)z, b(-L~) - u{-L~)z) - g 2 (a(-L _ ), b(-L~)) 
= e + (O(e) + zh v {z))z. 
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The functions h$ and h v are analytic on \z\ < (K + l)/i and we have h^(z) = O z (l) and 
h v (z) = O z (l). Substituting £ and rj above into (4.13a), we obtain 

X = (1 + 0(e) + zh^{z))z - ^(9; fi) - fiW u (e - ^(9; fi) + (0(e) + zh v (z))z, 9; /i) 

= (1 + 0(e) + zh^z))z - ii<f>(6; n) - fiW u (e - ^(9; fi), 9; n) 

- fiW u (e - ^{9- fj,) + (0(e) + zh v (z))z, 9; fx) + fiW u (e - ^(9; fi),9; fi). 

This implies 

(4.27) X = (1 + Oe^e) + zh(z, 9; n))z - ^(1) 

where h(z,9] /i) is analytic in z, 9, and \i and satisfies h = Zj e tfJ ,(l). Now substitute 

X = /iX, z = \iZ 

into (4.27) and note that |Z| < Kq + 1. We obtain the claimed formula for X. 
For the Y-component, we substitute £ and 77 above into (4.13b) to obtain 

Y = e + (0(e) + zh v (z))z - /-0 - ^W s ((l + 0(e) + zh^(z))z - /x0(0; //), 0; /x). 

Set Y = /jY and z = fiZ and note that \Z\ < K + 1. We obtain the claimed formula for 
Y. ■ 

Corollary 4.1. On E", we have 

Z = (1 + fljP (e) + /iO x ,, iP (l))(X + 9iP (l)). 

Proof. We start with (4.27). This equality is invertible and we have 

(4.28) z = (1 + Oe^e) + W/>(W, 0; //))W 
where 

W = X + KMl) 

and h(W,9;n) is analytic in W, 0, and /i and satisfies A = (9w,0,/z(l)- Writing (4.28) in 
terms of Z and X, we have 

Z = (1 + O e , p (e) + /iOx, e , P (l))(X + 0^(1)). 

■ 

Corollary 4.2. On we /iawe 

Y = /U - 1 £ + Xi , iP (l). 

Proof. We first regard Y as a function of Z, 0, and p using the formula for Y in Proposition 4.4 
and then regard Z as a function of X, 9, and p using Corollary 4.1. ■ 

Remark 4.3. Terms of the form /xOx,e, P (l) are not equivalent to terms of the form Ox,e, P (fJ>)- 
A term of the form nOx,o, p (l) has C 3 norm bounded above by K/i while a term of the form 
Ox,d, P (fJ-) has C 3 norm bounded above by K(e)/i. In estimates in Section 4.3B and 4.3C, we 
always have the former, not the latter. 

C. Conversion on E+. On S + we need to write X and Y in terms of Z. 
Proposition 4.5. On S + we have 

X = H- 1 e + O z ,eM 

Y = (1 + O e , p {e) + ^,, P (1))Z - O fljP (l). 



( 4 - 29 ) ,,r-s 
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Proof. On E+, s = L+. We have 

a(L + ) =e + Q 1 (e,0) = e + 0{e 2 ) 
b(L + ) = Q 2 (s,0) = O(e 2 
and 

(4.30) u(L + ) = -1 +0(e), v{L + )=0{e). 

Let (z, 6) G S + . We compute the values of X and Y for this point. Using (4.16) and (4.1), 
we have 

f = a (L + ) + v(L + )z + q l (a(L + ) + v(L + )z, b(L + ) - u(L + )z) 

= e + 0(e)z + qi (a(L + ) + v(L + )z, b(L + ) - u(L + )z) - qi (a(L + ), b(L + )) 

= e + (0(e) + zkz(z))z. 

Similarly, we have 

rj = b(L + ) - u(L + )z + q 2 (a(L + ) + v(L + )z, b(L + ) - u(L + )z) 

= -u(L + )z + q 2 (a(L + ) + v(L + )z, b(L + ) - u(L + )z) - q 2 (a(L + ), b(L+)) 

= {l + 0{e) + zk r) {z))z. 

We now write X and Y in terms of z using (4.13a) and (4.13b). The rest of the proof is 
similar to that of Proposition 4.4. ■ 

Corollary 4.3. If L + is sufficiently large, then Y > 1 on S + . 

Proof. This follows directly from the definition of S + . ■ 

5. Explicit computation of M and N 

In this section we explicitly compute the flow-induced maps M : E~ — > S + and N : S + — > 
E~. The map M : X~ — > E + is computed in Section 5.1. In Section 5.2 we study the time-t 
map of equation (4.15). The map N : E + — > E _ is computed in Section 5.3. 

5.1. Computing M : S _ — ► S + . Recall that s = —L~ on E~. Let q = (—L~,Z ,9 ) G 
XT and let (s(£), Z(t),9(t)) be the solution of system (4.22a)-(4.22c) initiated at the point 
Z , O )- Let i be the time such that s(t) = L + . By definition, M(q ) = (L + , Z(t), 6(t)). 
In this subsection we derive a specific form of M using (X, #)-coordinates to uniquely locate 
points on E~ and (Z, #)-coordinates to uniquely locate points on E + . Define 

K x (e) = - P A L eti +E ^ ds 

where 

L+ 

(u(s) + v(s))h(a(s),b(s))e- % E{r)dr ds 

' -L- 

is obtained by changing the integral bounds of the improper integral A in (2.5) to — L~ and 
L + . Also define 

p h — e f- L - E (s)ds_ 

Lemma 5.1. 

P L ~ < 1, K^e) ~ £~« > 1. 
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Proof. Both estimates follow directly from Lemma 4.2 and the fact that 



Proposition 5.1. Let (K ,9 ) £ £ and mie (Z, 0) = M(X ,# ). We /iai>e 
(5 ^ = #o + ^(L+ + L") + Ox ,eo,p(^) 

Z = Ki(e)(l + ci sin^o + c 2 cos# ) + Pl&o + ^W 1 ) + + XoA M) 

where c\ and c 2 are constants satisfying 



Proof. Using (4.22c), we have 

0(t) = 9 + cut. 

Integrating (4.22a) and (4.22b), for * e [-2L _ ,2L+] we have 

s (t) = —L~ + i + O t ,z ,e M- 
Inverting the last equality, we obtain 

t(s) = s + IT + O s ,z ,6 M- 
Substituting 9(t) and t(s) into (4.22b), we obtain 
dZ 

(5.2) — = E(s)Z - (u(s) + v(s))(pU(s, 0) + sin(£ + wL" + ws)) + a,z o ,0o,^). 

Note that in (5.2), (s, Z ,9 ,p) is such that s e [— 2L~,2L + ], (Z ,#o) £ and p = ln/i e 
(— oo,ln/i ]. Using (5.2), we obtain 

(5.3) Z{s) = P s ■ (Z - $ S (0 O ) + O s , Zo ,e M) 
where 

P s = 6 f-L- E ( T ) dT 

(5 4) r s 

$ s (9) = J (u(r)+v(r))(pm(T : 0)+sm(9 + ujL- + ujt)) ■ e~ E{f)df dr. 
From (5.3), it follows that 

9 = 9 + uj(L + + L-) + O Z(h eM 
Z = P L (Z -^ L+ (9 ) + ZoAh M). 

We want to write the right-hand side of (5.5) in (X , ^-coordinates. Using Corollary 4.1, 
we have 

9 = 9 + u;(L + + L-) + O M M 

Z = P L (X - $ l+ (0o) + O eo , p (l) + Ox O)0O) p(£) + XoAhP (fi)) ■ 
Let iT 2 be such that 

|Xq + OeU 1 ) + 0^e QyP {e) + Xo ,, OiP (/i)| < K 2 
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on and observe that by letting 

(5.7) K (e)= max 2 \P S (K 2 - $ a (0))| , 

0SS 1 
se[-2L",2L+] 

we conclude from (5.3) that all solutions of system (4.22a)-(4.22c) initiated inside of S~ will 
stay inside of 

{(s,Z,9) : s G [-2L-,2L + ], \Z\ < K (e)} 

before reaching s = L + . To finish the proof of Proposition 5.1, it now suffices for us to prove 
the following lemma. 

Lemma 5.2. For p e [pi,p 2 ], we have 

-P L $ L+ (6>) = #i(e)(l + ci sin 6> + c 2 cos 6>) 
where C\ and c 2 are constants satisfying 

Proof of Lemma 5.2. Recall that in (5.4), H(s, 0) = h(a(s), b(s)). We have 
P L $ L+ (6) = ef° + E ^ dT ■ ^ (u(s)+v(s))(ph(a(s),b(s)) 

+ sm{6 + uL- +us))e-f« E(T)dT ds 

= e# + e{t) dT ■ { P A L + (C L cos ujL- - S L sin ujL') sin 9 
+ (Sl cos ujL~ + Cl sin ujLT) cos 6) 



where 



A L = J {u{s) + v{s))h{a{s),b{s))e- ! ^ E(T)dT ds 

/L+ 
(u(s) + v(s)) cos(us)e~ ■!» E ^ dT ds 

,L+ 

Sl = J (u(s) + v(s))sm(us)e-fo E ^ dT ds. 



Observe that A, C, and S in (H2) are obtained by letting V s1 = oo in A L , Cl, and Sl- We 
now write 

(5.8) P L $ L + (9) = pA L eti + Eir) dT ■ (1 + ci sin 6> + c 2 cos 0) 



where 



(Cl cos uL — S l sin ujL ) 

Cl = 

(Sl cosu;L~ + Cl sinwL - ) 



c 2 



We have 



2 i 2 

Ci + c 2 



A L p 

(Cl + Sj) 
A\p> 
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Using (H2), for V s1 sufficiently large we have 

\A L -A\<±\A\ 



cj + s? - s/cpTs 



Therefore, for p G [pi,p2], where 

, , 202 VC 2 + S 2 396 VC 2 + 

(5 ' 9) Pl = -W^^' p2 = "ToT — a — ' 

we have 

\<4^ 2< \ 

Notice that because of the way in which p 1 and p 2 are defined, we have — pA L > 0. We also 
have 

K^e) = - P A L eti + E{T)dT ~ e~« 
from Lemma 5.1. Equation (5.8) for P L Q L +{9) is now in the asserted form. ■ 

By using Lemma 5.2, we can now rewrite (5.6) as (5.1). This finishes the proof of Propo- 
sition 5.1. ■ 

Remark 5.1. Observe that in formula (5.1) for Z, the term with Ki(e) in front dominates 
the second term because Ki(e) ^> Pl- The inclusion M(£ _ ) C S + follows directly from (5.1). 

5.2. On the time-t map of equation (4.15). The computation of N : S + — > £~ contains 
two major steps. The first step is to compute the time-t map of equation (4.15) inside 
U £ . This is done in Section 5.2. The second step is to compute the time it takes for a 
solution of equation (4.15) initiated in S + to reach E~. This is done in Section 5.3. These 
computations are technically involved because we need to control the C 3 norms of the map 
N on S + x (— oo, In p ], where the interval in the product is the domain of the parameter p. 

We start with the first step. Let VF(E + ) be a small open neighborhood surrounding S + 
in the space (X, Y, 6). In this subsection we let (X , Y , 9 ) G W(T I + ) and regard p = ln/i e 
(— oo, In po] as the parameter of equation (4.15). We study the time-t map of equation (4.15) 
assuming that up to time t, all solutions initiated from W{Y> + ) are completely contained 
inside U £ . Recall that in equation (4.15), F(X, Y, 9; p) and G(X, Y, 9; p) are analytic on 

D = {(X, Y, 9, p) : p e [0, no], (X, Y, 9) e UJ 

where 

U e = {(X,Y,6): ||(X,Y)|| <2ep-\ 9 G S 1 }. 
Forg = (X ,Y ,^o)G^(S + ),let 

q(t,q ;p) = (X(t,q Q ; p),Y(t,q ; p),9(t,q ; p)) 
be the solution of equation (4.15) initiated from q at t — 0. Using (4.15), we have 

X(t, q ; p) = x e^ { - a+ ^ {q(s ' qo ^ M)ds 
(5.10) Y(t, q ; p) = Y e^ i/3+fMGiq{s ' qtr '^ ))ds 

9(t, q ; p) = 9 + cot. 
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We now introduce the functions U(t, q ; /x) and V(t, qo', A*) and rewrite (5.10) as 

X(t, q ; fi) = X e^ a+u ^ qo '^ t 
(5.11) Y(t, g ; AO = Y e (/3+v( *' 9o;/i)) * 

6>(t, g ; Ai) = 6»o + art. 



Using (5.11), we have 
(5.12) 

We also have 
(5.13) 



U(t,q ;n) = t In + a 

V(t,q ;n) = t In /3. 



U{t,q ; fi) = t 1 / /j,¥(q(s,q ;n);ii)ds 
Jo 

V(t,q ;n) = t~ 1 / fiG(q(s,q ;fi);fi)ds. 
Jo 

In the next proposition we regard U = U(t,q ; fx) and V = V(t,qo; fi) as functions of t, qo, 
and p and we write U = U tAOtP and V = V tMhP , respectively. We define the domain of these 
two functions as follows. Let 

mW = {Qo e W(Z + ), p G (-oo,ln/i ], t G [l,T(q ,p)}} 

where the upper bound T(q ,p) on t is designed to keep the solution inside U e . 

Proposition 5.2. TTiere exists K > suc/i £/iat 

ll#t,</o,pllc3(B t ,, , p ) < ll^,</o,pllc3(D t ,, , p ) < KfX. 

Proposition 5.2 is proved in Section 7.1. 

Remark 5.2. By combining Proposition 5.2 and (5.11), we can now write the time-t map 
from W(E+) to U £ as 

X(t, Xq, Y , 6» ; /x) = Xoe^- a+0 ^oWo, P M)t 
(5.14) Y(t, Xo, Y , 9 ; fx) = ¥ e ( ^°'^^- 9 o^' l))i 

0(*,Xo,Y o ,0o;aO = O + ojt. 

5.3. Estimates on T(Z ,9o,p). For g = (-Zo,#o) G S + , let q(t,qo',fi) be the solution of 
equation (4.15) initiated at g and let T be the time this solution reaches S~. In this 
subsection we regard T as a function of Zo, 9q, and p and we obtain a well-controlled 
formula for T that is explicit in the variables Zo, 9q, and p. Since the images of M are 
expressed in (Z, #)-coordinates through (5.1), we must write the initial conditions for Ji in 
(Z, ^-coordinates on £+ to facilitate the intended composition of N and M. 

Estimates on T(Z ,9o,p) are complicated partly because as a function of Z and 9q, h 
is implicitly defined through equations written in (X, Y, #)-coordinates on T, ± . The com- 
putational process therefore must involve (5.14) and the coordinate transformations on S ± 
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presented in Sections 4.3B and 4.3C. Before presenting the desired quantitative estimates, 
we explain how to obtain T(Z ,9 ,p) in a conceptual way. Using (5.11), we obtain 

X(T, Xq, Y , 9 ; fi) = X e^ a+U ^ M ^ T 
(5.15) Y(T, Xo, Y , 9 ; f*) = Y e^ +V( ™^°^ T 

9(T,X ,Y ,6 ;ii) = e + uT. 

In (5.15), X and Y are not independent variables. These quantities satisfy 

Yo = (1 + O d(hP (e) + /i0 ZO)flbiP (l))Z o - C\, p (l) 
by Proposition 4.5. We write 

X(T)=X(T,Xo,Y Q ,9 Q ;n) 
Y(T)=Y(T,X ,Y ,e ; f i) 
6{T) = 6 + ujT. 

By definition, X(T), Y(T), and 0(T) are also related through Corollary 4.2. For the benefit 
of a clear exposition, we write the conclusion of Corollary 4.2 as 

Y = eyT x + f(X, 9;p) 

where 

f(X,0;p) = 0x^(1). 

We have 

(5.17) Y(T) = £/ x- 1 + f(X(T)^(T);p). 
We use (5.15) to implicitly define T(Z ,9 ;p). We have 

(5.18) Y(T) = Yoe (/3+v(T,x ,¥ , eo;P ))T > 

The right-hand side of (5.18) is relatively simple: we only need to substitute for X and 
Y using (5.16). The left-hand side of (5.18) is conceptually more complicated. We need to 

(1) Write Y(T) as a function of X(T), 9(T), and p using (5.17). 

(2) Substitute for X(T) and 6*(T) using (5.15), thereby obtaining Y(T) in terms of T, X , 
Y , ^o, and p. 

(3) Use (5.16) to write X and Y in terms of Z and 9 . 

After all of these substitutions are made, we regard (5.18) as the equation that implicitly 
defines T(Z , 9 ;p). We use this equation as the basis for the computation of T(Z , 9 ;p). 

Proposition 5.3. As a function of Z , 9 , andp, the map T satisfies 

HT-^ln^- 1 !!^ < K. 
Proposition 5.3 is proved in Section 7.2. 
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5.4. Computing N : £+ — ► £ _ . We derive a formula for the induced map 3\f p : S + — > E _ . 
For (Z , 6> ) G £+, we write (X l5 ^) = K P (Z , 6> ). We start with [/ and V in (5.15). 

Lemma 5.3. On S + x (— oo, ln// ], we /jat>e 

U(T,X ,Y ,9 ;p) = fJ-OzofioAV 
V(T,X ,Y ,e ;p) = nOz ,8 ,p(l)- 

Proof. We write J7 and V as functions of (Zq,9q,p) using Proposition 5.3 for T(Zq,9q;p) 
and (5.16) for Xo and Yo- This lemma is established by applying the chain rule and using 
Proposition 5.2, Proposition 5.3, and (5.16). ■ 

Proposition 5.4. The flow-induced map N p : S + — > S~ zs gwen 6?/ 



(5.19) \£ + ^OzofioA 1 ) 



-i 



X i = 7TT (t 1 + O<w00 + ^ Zo ,, 0iP (l)]Z - O^l))^ 



= O + u l n . ( £ + ^o,e , P (l))^ 1 



P + »O Z0fi0tP (l) [1 + Og 0tP (e) + ^0 ZoAhP {l)]Z, - 9(hP (l) 



where 



a = a + nO Zofio , P { l )i P = P + ^OzofioA 1 )- 
Proof. Using (5.17), (5.18) and Lemma 5.3, we have 

1 ln Y(T) 



,, 2m /3 + ^Oz oA , p (l) F 

1 • ] i (e + zifcx^^cr)^))//- 1 

/? + ^O Z0i , 0iP (l) 11 Y 
By using Proposition 5.3 and the fact that f(X,9;p) = 0x,0,p(l), we have 

f(X(T)^(T);p) = O^ 0)P (l). 

Now (5.20) gives 

( 521 ) T _ 1 ln /i" 1 (£ + ^C Z()A)iP (l)) 



/5 + vO ZoAhP (l) [1 + a ,p(e) + fiO ZoAhP (l)]Z - O eo M 

Here we use (5.16) for Y . 

The desired formula for 9 1 now follows from 9\ = 9 + uoT. For X x we use 

Xl = fi-\e + fiO ZoA) M)e- {a+ ^oA^ 
and substitute for T using (5.21). ■ 

6. Proof of Theorem 2.1 

In Subsection 6.1 we compute 5F p = N o M by using Propositions 5.4 and 5.1. In Sub- 
section 6.2 we apply the theory of rank one maps to the family {9^}, thereby proving the 
existence of rank one chaos as claimed in Theorem 2.1. 
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6.1. The flow-induced map = No M. We regard p as the fundamental parameter of 
the flow-induced map 5 : E - -> E". For (X O ,0 O ) e let (X^) = pSf o M)(X , 9 ). We 
compute SFp : (X ,#o) ^ (Xi,#i) by combining (5.19) and (5.1). 

Proposition 6.1. The map H p : E~ — > E~ groen fry 



(6.1a) X 1 = (/i(£ + 0XoA ) ,p(/^))" 1 )' 3 

x ((! + + OxoA»)Z " <£(!)) ' 

(6.1b) 9 l = 9 + u;(L+ + L") + Xo ,0 O)P (/i) 



+ Xo , eo M (1 + Ogj^(e) + O m ,M)Z - 
w/iere 

Z = #i(e)(l + ci sin^ + c 2 cos# ) + ^l[X + O 6o , v {l) + C Xo ,eo,P^) + OxofioM] 

a = a + Xo ,e , P (^) 
(3 = (3 + O Xafi0 M 

and the superscript <9X (/i) on a given term indicates that the partial derivative of the term 
with respect to X is C(/i). We also have 



tfi(e) ~e «, - < ^c? + cl < -. 

Proof. We first examine the formulas for a and (3. The error terms in Proposition 5.4 have 
the form 

and Z and 9 are given in terms of X , 6> , and p by (5.1). Using (5.1), we see that the C 3 
norms of Z and 9 are < K{e). It follows from the chain rule that 

We follow the same line of reasoning to compute X x and 9\. We replace Z and 9q with Z 
and 9 in (5.19) and then substitute for Z and 9 using (5.1). Using (5.19), we have 



( , 2) ^M^w' (|i + ^) + ^,i)„- OU1) ) 



0! = + 7; ^ 7TT In 



+ ^W 1 ) [1 + O^(e) + ^,(1)]^ - 0^(1) ' 

In (6.2), terms of the form ^O^eJS) are rewritten in the form 0x o ,0o,p(aO using (5.1). 
Terms of the form 0§ (e) are rewritten in the form 0^g^ p (e) because the C 3 norm of 9 is 
bounded by a constant K independent of e and because = 0([i). Reasoning analogously, 
terms of the form 0§ (1) are rewritten in the form O^ g^ p (l). ■ 
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6.2. Proof of Theorem 2.1. We are finally ready to prove Theorem 2.1. 

The two-parameter family {3a,&„}- We write {3^} as a two-parameter family { of 
2D maps. Both a and b n are derived from \i = e p as follows. Let /io > be sufficiently small. 
Define 7 : (0,/io] - ► K via j(fj) = ^ln/i _1 . For n G Z + satisfying n ^ (27r/3) _1 u; In //q x , let 
(j, n G (0, /io] be such that 7(/x n ) — n - Notice that //„ — > monotonically. Set 6 n = /i n . For 
// G (fJ-n+i, (J>n\ and a G [0, 27r) = S 1 , we define 

fi(n, a) = 7~ 1 (7(/i„) + a) = ji n e ^ a 



and 



o) = In n(n, a) — In fj, n — —a. 



Define 

$a,b n — 3p(n,a)- 

Verification of (CI)— (C4). We prove Theorem 2.1 by applying Propositions 3.2 and 3.3. 
We verify (C1)-(C4) for ^ a ,b n - Proposition 6.2 establishes (CI). 

Proposition 6.2. We have 

(6.3) ||9 r a, 6 „(X^)-(0,^ (X^))|| C 3 (s - x[0 ^ )) ^0 
as b n — > 0, where 

5F a ,o(X, 0) = + + L") + a + ^ lnfe^e)" 1 ) 

(6.4) - - In 



P 



1 + e , p (e)) ^1 + ci sin(0) + c 2 cos(#) 

+ (X + a, P (l) + ^(e)) ) - K^OoM 

Proof. The only problematic term in (6.1b) has the form 

^ ln -1 

which we write as 

Observe that the C 3 norm of the second term — > as & n — > and the first term may be 
computed modulo 2-7T and is therefore equal to a. Viewing /iasa function of a, the C 3 norm 
of Xi is bounded by 

and therefore decays to as b n — > provided that (Hl)(b) holds. ■ 
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For (C2) we apply Proposition 3.1 to the family of circle maps 
^0,0(0, 9) = 9 + lu(L+ + L-)+a+^ ln^^e)" 1 ) 



(6.5) 



^ln 



(1 + O e , P {e)) [I + ci sin(0) + c 2 cos(6>) 



To apply Proposition 3.1 to the family {3^0(0, #)}, we set 

V(9) = - ln(l + ci sin 6» + c 2 cos 0) 
$(0, a) = 5F a , (0, 0) - 7 - - a - XV(9) 

where 

7 = u;(L + + L-) + ^ln( £ K 1 ( £ )- 1 ). 

The assumption on the C 3 norm of $ is satisfied if e is sufficiently small. 

Hypothesis (C3) follows directly from (6.4). Hypothesis (C4) follows from a direct com- 
putation using (6.2). Finally, to apply Proposition 3.3 we need to verify that Ao > In 10. 
This follows if ui is sufficiently large. The proof of Theorem 2.1 is complete. 

7. Computational proofs 

7.1. Proof of Proposition 5.2. Let F = F(X, Y, 9; //) and G = G(X, Y, 9; /x) be as in 
equation (4.15). For a combination Z = ~K dl Y d2 n d3 of powers of the variables X, Y, and /i, 
let d\ denote the corresponding partial derivative operator, where k = d\ + di + d 3 is the 
order. There exists K 3 > such that for every Z of order ^ 3 and ^ i ^ 3, we have 

(7.1) |0*(^F-Z)| <K 3 , |d*(^G-Z)| <K 3 

on U)t,g ,p. This is because the C 3 norms of F(X, Y, 9; ji) and G(X, Y, 9; /j) are bounded on 
U £ x [0,'no] and because F(X, Y, 9; /j) = F(//X, /iY, 0; fi) and G(X, Y, 0; fj) = G(/iX, /iY, 0; fi). 
C° estimates. Using (7.1) with i = k = and (5.13), we have 

(7-2) II^I|co(b Mo , p ) < II^Hcocd^,,) < K 3 fl. 

C 1 estimates. We now estimate the first derivatives. 

On dY U and <9y U. Using 9(t) = 9 + wi, we have 8y 9 = 0. Using (5.13), we have 
(7.3a) dY U = /it" 1 / {dxF • <9y X + ^F • 8y Y) rfs 

(7.3b) d Yo V = nr 1 [ (dxG-dY X + dYG-dy Y)ds. 

Jo 

To make these formulas useful, we need to write <9y X and <9y Y in terms of 8y U and <9y U. 
For this purpose we use (5.12). We have 

d Yo X = tXd Yo U 

(7.4) Y 

d Yo Y = tYd ¥o V + —. 

-U-0 



34 QIUDONG WANG AND WILLIAM OTT 

Combining (7.3a), (7.3b), and (7.4), we obtain 



(7.5) 



d ¥o U = fit' 1 [ (<9xF • X • sdY U + (9yF • Y • s^V) ds + fit' 1 [ ^Y-^-ds 
Jo Jo ¥ o 

dy V = ^r 1 / (9xG • X • sdY U + <9yG ■ Y • sd^F) + fir 1 [ <9yG • — ds. 

Jo Jo ^o 



Using (7.1), we have 

X\<K 3 , \dxG-X\ < K 3 , |dyF-Y|<Ar 3 , \dyG ■ Y| < K 3 . 



Using (7.5), we have 



(7.6) 



\dy U\ < Kfrt- 1 f (\sdY U\ + \sdY V\) ds + Kfi 
Jo 

\dY V\ ^ Kfir 1 [ (\sdY U\ + Isdy.Vl) ds + Kfi, 
Jo 



from which it follows that 

\dy U\ < Kfi, \d* Q V\ < Kfi. 

On dx U and dx V. Mimic the proof above. 

On dg U and dg a V. We follow similar lines of computation. Since dg 9 = 1, we have 

d 0o U = fir 1 [ {dyj ■ d 0o X + SyF • d do Y + d e ¥) ds 
Jo 

d 9o V = fir 1 [ (dxG ■ d 0o X + <9yG • dg Y + d e G) ds. 
Jo 

Analogous to (7.4), we have 

d eo X = tXd 6o U, d 9o Y = tYd 6o V. 
Arguing as above, we conclude that 

\d eo U\ < K/i, \d eo V\ < Kfi. 

On d p U and d p V. We follow similar lines of computation. Note that we have 

d p fi = fi, d p W = fidpW, 
and so on. Starting with (5.13), we have 



(7.7) 



d p U = fit' 1 [ ¥ds + fir 1 [ (<9 X F • d p X + 8y¥ ■ d p Y + fidpF) ds 
Jo Jo 

dpV = fir 1 [ Gds + fir 1 [ (5xG • d p X + 8yG ■ d p Y + fid^G) ds 
Jo Jo 

and using (5.12) we have 
(7.8) 



d p X = tXdpU 
d p Y = tYd p V. 



Now argue as above. 

On d t U and d t V. The partial derivatives of U and V with respect to t are easier to estimate 
because when differentiating with respect to t using (5.13), no derivatives are involved on 
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the right-hand side so the estimates on d t U and d t V are obtained directly from C° estimates. 
We have 

\d t U\ < Kfi, \d t V\ < Kfi. 
This completes the desired estimates on the first derivatives. 

C 2 estimates. We now move to the second derivatives. We estimate d YoYo U and <9 2 „ V „U 
first. Using (7.3a), we have 

9 2 y y U = fir 1 f (c^F • (d^X) 2 + 2dl Y ¥ • (9 Yo X)(a ¥o Y) + ^(F • ^Y) 2 ) ds 
Jo 

+ fit- 1 [ (d x ¥ ■ a 2 o¥o X + <9yF • 9 2 o¥() Y) ds. 
Jo 

Using (7.4), we have 

% oYo x = td yo x-d Vo u + tx% oYo u 

d 2 YoY J = tck.Y -^V + tY- di oYo V + %^ - 
Therefore, d Y Y U is given by 



(7.10) 



$y y U = M*" 1 f ■ (9y ^) 2 + 2di ¥ F • (^ X)(av Y) + dw(F • d^Y) 2 ) ds 

Jo 

+ fit' 1 J (3xF • d* X ■ sdy U + dyF • dy Y • sdy V) ds 

+ fit' 1 [ (&F • X • sdl )Yo U + dyF • Y • sdl )Yo V) ds. 
Jo 



To estimate the first three integrals in (7.10), we use (7.4) for <9y X and 9y Y. Using the first 
derivative estimates and using (7.1) repeatedly, we bound these integrals by Kfi. Note that 
we also need Yo > 1 (see Corollary 4.3) for the third integral. Together with an analogous 
formula for d\ Y V in which we replace F with G, we conclude that 



dl Y U\ < Kfi, \d 2 V\ < Kfi. 



All other second derivatives are estimated similarly. Here we skip the details to avoid 
repetitive computations. 

C 3 estimates. Third derivatives are estimated in the same spirit. Since the formulas for a 
given third derivative depend on previous computations of relevant second derivatives, here 
we estimate dy oYop U and d YoYop V as a representative example. Of all of the third derivatives, 
these are the most tedious to compute. 

To compute 9 YqY U we apply d p to (7.10). The explicit factor fi written in front of all 
integrals generates a collection of terms that is identical to the right-hand side of (7.10). We 
showed when estimating second derivatives that the size of each of these terms in bounded 
by Kfi. 

The remaining terms are produced by applying d p to the functions inside of the integrals 
in (7.10). The terms produced from the first three integrals are estimated using the C 2 
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estimates. Estimate (7.1) is used repeatedly. It is critically important that potentially prob- 
lematic terms in the form of powers of Y and X, introduced by using the likes of (7.4), (7.8), 
and (7.9), are always matched perfectly with corresponding partial derivatives with respect 
to F or G. Applying d p to the fourth integral, we obtain an integral term of the form 

(I) = lit- 1 f {djJ-X- sd YoYop U + c\¥ • Y • sd* oYop V) ds 
Jo 

and a collection of other terms that can be treated the same way as the terms produced by 
differentiating the first three integrals. We have 

< Kf,r l f (K oYop C/| + \sdl Yop V\) ds. 
Jo 

Combining this analysis with analogous estimates for |<9 Yo¥op U|, we obtain 

,U\<KfM, \dl YnP V\<Kfi. 



This completes the proof of Proposition 5.2. 

7.2. Proof of Proposition 5.3. The proof of this proposition is lengthy because of the 
complicated composition process explained earlier in Sect. 5.3. 

C° estimates. We first establish a C° control on T. 

Lemma 7.1. There exist constants K4 < K 5 independent of e such that for all q e S + , we 
have K4 In [T 1 < T(q ; /i) < K 5 In . 

Proof of Lemma 7.1. Using 

Y(T) = Y e^ +V ^ T 

we obtain 

T — \ In^gl. 
P + V(T) Y 

Since (X(T), Y(T), 0(T)) is on E~, Proposition 4.4 implies that 

Y(T) « // _1 £ 

and the desired estimates follow from | V(T") | < Kfi and 1 < Y < K(e). ■ 
Lemma 7.2. We /jai>e /j,~~ 1 e~ aT < 1. 
Proof of Lemma 7.2. We substitute 



T = \ In^ 

/3 + V(T) Y 



into (5.15) to obtain 



X(T) 



q-t/(T) 

Y \ "+ V < T ) 



y(t) ; 

We then use Y(T) « s^ 1 , X « e/x -1 , |f/(T)| < K/i, \V(T)\ < Kfi, and a > f3 to conclude 
that X(T) < e. We have 

^ -1 e- QT < X e(- a+l/ ( T )) T = X(T) « e. 

For the first inequality, we use X rs e/x -1 and |[/(T)T| < Kfilnfi" 1 <C 1. This proves the 
lemma. ■ 
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C 1 estimates. We present C l estimates with respect to (Z ,9 ,p), where (Z ,9 ) £ S + 
and p G (—00, In fi ]. 

Lemma 7.3. There exist constants K 7 and K 8 independent of e such that 

||X(T)|| c i <K 7 + K 8 \\T\\ c i, \\e(T)\\ c i<K 7 + K 8 \\T\\ c i. 

Proof of Lemma 7.3. The bound on 9{T) is trivial because 9{T) = 9 + ojT. For X(T), we 
have 

X(T) = X e^ a+U ^ T 

= £/[x - 1 e (- a + l/ ( r . x o,¥ ,e ;p))T + ^ e (-a+f7(T,X o ,Yo,0 o ;p))T_ 

Notice that for the second equality, (5.16) is used for X . We regard X and Y as functions 
of Z , 9 , and p defined by (5.16). The desired estimate follows from using Proposition 5.2 
for U and (5.16) for X and Y . We also use Lemma 7.2. ■ 



Lemma 7.4. We have 



|T--ln^- 1 || c i < K. 



Proof of Lemma 7.4- Using (5.17), we write (5.18) as 

p-\e + /if (X(T), 9(T); p)) = Y e^ +V ^ T . 

Solving for T, we obtain 

( , n) ^^-w^w^-j^vm^ 

+ ,3 + vtr) ln(£ + ' ,f (x(r) ' eiT ^- 

In (7.11), V(T) = V(T,X ,Y ,9 ;p), and X and Y are written in terms of Z , 9 , and p 
using (5.16). Using Proposition 5.2, we have 

HT-iln^Hco < K. 

First derivatives of T are estimated by directly differentiating (7.11). We estimate dz T 
as a representative example. Differentiating (7.11), we have 

d Zo T = (I) + (II)d Zo T, 

where (I) is a collection of terms that do not depend on d Zo T and (II) is a function of Z , 9 , 
and p. Using Proposition 5.2 for V(T), (5.16) for X and Y , and Lemma 7.3 for c^ X(T) 
and d Zo 9(T), we have < K and \{II)\ < 1. 1 ■ 

Higher derivative estimates. With the first derivatives controlled by Lemmas 7.3 and 7.4, 
we estimate the second derivatives by first proving a version of Lemma 7.3 and then proving 
a version of Lemma 7.4 for the C 2 norms. We then do the same for the C 3 norms. This 
completes the proof of Proposition 5.3. 
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